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A separable form of pairing interaction in the channel has been introduced and successfully 
applied in the description of both static and dynamic properties of superfluid nuclei. By adjusting 
the parameters to reproduce the pairing properties of the Gogny force in nuclear matter, this sepa- 
rable pairing force is successful in depicting the pairing properties of ground states and vibrational 
excitations of spherical nuclei on almost the same footing as the original Gogny force. In this ar- 
ticle, we extend these investigations for Relativistic Hartree Bogoliubov theory in deformed nuclei 
with axial symmetry (RHBZ) using the same separable pairing interaction. In order to preserve 
translational invariance we construct one- and two-dimensional Talmi-Moshinsky brackets for the 
cylindrical harmonic oscillator basis. We show that the matrix elements of this force can then be 
expanded in a series of separable terms. The convergence of this expansion is investigated for various 
deformations. We observe a relatively fast convergence. This allows for a considerable reduction in 
computing time as compared to RHBZ-calculations with the full Gogny force in the pairing channel. 
As an example we solve the RHBZ equations with this separable pairing force for the ground states 
of the chain of Sm-isotopes. Good agreement with the experimental data as well as with other 
theoretical results is achieved. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.60.-n, 21.10.Dr 



I. INTRODUCTION 



Covariant density functional theory (DFT) based on 
Relativistic Mean Field (RMF) theory provides a mi- 
croscopically consistent description of the nuclear many- 
body problem l]. Conventional DFT without particle- 
particle (pp) correlations can be applied only for a few 
doubly closed shell nuclei. For the vast majority of nu- 
clei and in particular those far away from the /3-stability 
line, that play an important role in astrophysical appli- 
cations, the inclusion of particle-particle correlations is 
essential for a quantitative description of many phenom- 
ena in nuclear structure. In the framework of DFT pair- 
ing correlations are taken into account in the form of 
Hartree-Bogoliubov theory 0, d, 0, [E 111 for the ground 
states. Because of the numerical complexity of the de- 
formed Hartree-Bogoliubov equations, monopole pairing 
or density dependent (5-pairing interactions have been 
widely used in the literature for deformed mean field cal- 
culations [5,]. But the results are affected by a cutoff 
parameter which has to be introduced in a rather arbi- 
trary way. Recent investigations show clearly that many 
results depend on this cut-off Q- Therefore, in order 
to avoid the complicated problem of a pairing cutoff, the 
finite range Gogny force |8[ has been applied in many rela- 
tivistic applications 0] . The parameters of this force have 
been adjusted very carefully in a semi-phenomenological 
way to characteristic properties of the micro scop ic effec- 
tive interactions and to experimental data jlOl . [T]| . It 
has been shown in several applications that the RHB 
model with Gogny pairing provides an excellent tool for 
the description of ground state properties in finite nuclei. 
Such investigations have been first devoted to spherical 



applications, such as to halo phenomenon in light nu- 
clei [isl. to the properties of nuclei near the neutron 
drip line [14| , to the reduction of the spin-orbit potential 
in nuclei with extreme isospin values , to ground-state 
properties of Ni and Sn isotopes [16[, etc. RHB theory 
with Gogny pairing has also been applied to investigate 
deformed [itI, [H [li, [l^] and rotating 20] nuclei, but 
such calculations are limited due to their numerical com- 
plicity. 

Recently we have introduced a new separable form of 
the pairing force for RHB calculations in spherical nu- 
clei [2l|, [22| ■ The parameters of this separable force are 
adjusted to reproduce the pairing properties of the Gogny 
force in nuclear matter. It preserves translational invari- 
ance and has finite range. A similar ansatz has been 
used in the pairing channel of non-relativistic Skyrme 
calculations in Refs. jH, H^l- This pairing interaction is 
separable in momentum space. In r-space the transla- 
tional invariance leads to a (5-force in the center of mass 
coordinates and therefore, at a first glance, translational 
invariance forbids exact separability. However using well 
known techniques of Talmi and Moshinsky [2^, [2^, [23] it 
has been shown in Ref. [2lj that this force can be rep- 
resented by a sum of separable terms which converges 
quickly. This avoids the complicated problem of a cut- 
off at large momenta or energies inherent in zero range 
pairing forces. As we discussed in Fig. 6 of Ref. [2l|, we 
found that although the (5-force can give the same average 
gap as the Gogny force DIS if the size of the strength is 
adjusted properly, the individual matrix elements of the 
forces and the matrix elements of the pairing field A are 
very different from each other. The (5-force behaves in 
some sense very much like a constant pairing. In con- 
tract our separable force has a very similar behavior to 
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the Gogny force. 



II. THEORETICAL FORMALISM 



This simple separable force can reproduce the pairing 
properties of the ground-state for spherical nuclei on al- 
most the same footing as the original Gogny pairing in- 
teraction. Recently it has also been applied for studying 
dynamic properties of spherical nuclei [28| . The relativis- 
tic quasiparticle random phase approximation (RQRPA) 
based on the same separable pairing interaction in the 
pairing channel was used for calculations of low-lying 2+ 
and 3~ excited states in a chain of Sn-isotopes, which 
are very sensitive to the pairing channel. In comparison 
with experimental data and with the results of the origi- 
nal Gogny force it was shown that this simple separable 
pairing interaction is also very successful in depicting the 
dynamical pairing properties of vibrational excitations. 



So far, the new separable pairing force has been used 
only in spherical nuclei. In this article we apply this 
separable pairing force for axially deformed relativistic 
Hartree-Bogoliubov (RHBZ) calculations. For axially 
symmetric shapes the densities are invariant with respect 
to a rotation around the symmetry axis, which is taken 
to be the z-axis. Therefore it is convenient to work in 
cylindrical coordinates. The Talmi and Moshinsky tech- 
niques [1^ m, ll^] used in Ref. [2lj are restricted to 
spherical coordinates. Therefore we had to develop sim- 
ilar techniques for cylindrical coordinates working in an 
anisotropic oscillator basis. Again the matrix elements 
of this pairing force in this basis are no longer fully sep- 
arable. However they can be expanded, as in the spher- 
ical case, in a series of separable terms. Obviously the 
convergence of this expansion is not as fast as in the 
spherical case, it is still quick enough to save consid- 
erable numerical effort as compared to the full Gogny 
calculations. Finally we investigate for a chain of Sm- 
isotopes the ground-state properties such as binding en- 
ergies and deformations using the RHBZ-program with 
this new separable pairing force. As it is known, the 
ground-state properties of the deformed nuclei are highly 
affected by the pairing gap. Good agreement is found, 
when comparing with experimental data and with the- 
oretical calculations using RMF + BCS theory with a 
constant pairing force, or using non-relativistic Hartree- 
Fock Bogoliubov (HFB) theory with the full Gogny force 
DIS. 



We start our investigations in symmetric nuclear mat- 
ter with various densities. The gap equation in the ^Sq 
channel has the form, 

where a separable form of the pairing force is intro- 
duced fH, 



{k\Kfp"\k') = -Gp{k)p{k') 



(2) 



A simple Gaussian ansatz p{k) = is assumed. In 

Ref. [21] the two parameters G and a have been fitted to 
the density dependence of the gap at the Fermi surface 
A(fci?). Comparing with the Gogny force, we found two 
sets of parameters G = 738 MeV-fm"^ and a — 0.636 fm 
for the parameter set Dl [1] and G = 728 MeV-fm^ and 
a = 0.644 fm for the set DIS 

In the Hartree approximation for a consistent mean 
field, the RHB equations [1, |3| read 



ho-X A \ /C/,(r)\ fUkir) 



where h]j is the single nucleon Dirac Hamiltonian, 



(3) 



hD=~ia- V -I- (3{m + ga(j{r)) + g^TsUJ^ir) + gpP°{r) 

(4) 



The Dirac Hamiltonian contains the mean-field poten- 
tials of the isoscalar scalar a-meson, the isoscalar vector 
tj-meson, the isovector vector p-meson, as well as the pho- 
ton, m is the nucleon mass and the term —m subtracts 
the rest-mass and normalizes the energy scale to the con- 
tinuum limit. The chemical potential is to be determined 
by the subsidiary particle number condition, where the 
expectation value of the particle number operator in the 
ground state equals the number of nucleons. The col- 
umn vectors are the quasiparticle spinors and Ek are the 
quasiparticle energies. A is the pairing fields, which is 
an integral operator with the kernel 

Aabir, r')^^Yl ^abcdir, r')K,d{r, r') , (5) 



The paper is arranged as follows. The theoretical for- 
malism of RHBZ with the separable form of the pairing 
interaction is presented in Sec. II. The convergence of the 
expansion of the pairing force elements in the cylindri- 
cal harmonic oscillation basis is investigated in Sec. HI. 
In Sec. IV the ground-state properties of a chain of Sm- 
isotopes are calculated in the RHBZ approach. They are 
discussed in Sec. IV. Finally we give a brief summary in 
Sec. V. 



where a, 6, c, d denote the quantum numbers that specify 
the Dirac indices of the spinor. They run over the two 
spin orientations and the large and small components. 
Vabcdi^i^') are matrix elements of two-body pairing in- 
teraction. In general this should be a relativistic force 3 
and involve large and small components. However, since 
pairing correlations in nuclei are a purely non-relativistic 
effect it has been shown in Ref. [2^ that we can ne- 
glect the pairing matrix elements between large and small 
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components as well as the effect of the pairing matrix el- 
ements between small components and we consider only 
the upper part of the pairing field A in Eq. ([3|). The 
pairing tensor is defined as 



K,d(r,r') = V UckiryVdkir') 



Ek>0 



(6) 



For the axially symmetric deformed shape rotational 
symmetry is broken and therefore the total angular mo- 
mentum J is no longer a good quantum number. How- 
ever, the densities are still invariant with respect to a 
rotation around the symmetry axis, which is taken to be 
the z-axis. It then turns out to be useful to work with 
the cylindrical coordinates 



(7'^ cos (fi, rj_ sin (p, z) 



(7) 



In these coordinates the Dirac equation can be reduced 
to a coupled set of partial differential equations in the 
two variables z and r±_ that are solved by an expansion 
in an anisotropic harmonic oscillator basis [30l |. 

Since the interaction in the particle-hole (ph)-channel 
is identical to earlier calculations, here we only discuss 
the derivation of the matrix elements of the pairing inter- 
action for the separable form of Eq. Q in the pp-channel. 
First, we transform the separable force in Eq. Q from 
momentum space to coordinate space and obtain 



^('"i,'"2,r3,r-4) = 



-G5{R-R')\{l-P.)P{r)P{r') 



(8) 



where R ~ ■^('"i + '"2) and r = ri — r2 are the center of 
mass and relative coordinates of two paired particles, re- 
spectively. And P{r) is obtained from the Fourier trans- 
form of p(fc). 



P[r) = P{z,r±_) 



(47ra2)3/2 



(9) 



The term 5{R — R') in Eq. ^ insures the translational 
invariance. It also shows that this force is not completely 
separable in coordinate space. 

We start from the eigenfunctions of the deformed har- 
monic oscillator: 

|a) = \nz,nr,mi,ms) 

^yb;,b\ bz bi_ V27r 

with 



(11) 



The quantity b^ = ^/h/rmoz and b±^ — y^h/muj± are the 
harmonic oscillator length. The polynomials Hn{C) and 
L™{ri) are Hermite polynomials and associated Laguerre 



polynomials as defined in Ref. [3l|. The normalization 
constants are given by 



1 



[ur + |mi|)! 



(12) 



where the m; and vrig are the components of the orbital 
angular momentum and of the spin along the symmetry 
axis. The eigenvalue of j^, which is a conserved quan- 
tity in these calculations, is = to; + m^, the parity is 
given by TT = (— )"^+'"!. The conventional deformation 
parameter (3 is obtained from the calculated quadrupole 
moments through 



Q = Qn + Qp = \j^-^^ARl[3 , (13) 

with _Ro = 1.2 A^/'^ fm. Where The quadrupole Qn,p 
moments for neutrons and protons are calculated using 
the expressions 

Qn,p = (2r2p2(cOS0))„,p = (2z2 - - l/)n,p . (14) 

In the pairing channel, although the total angular mo- 
mentum J is no longer a good quantum number, we still 
have ill + f72 = 0. Together with the projector ^{1 — Pa), 
the two-particle wave function can be written as: 

|ai,Q;2) = |n^i,n^2)|nriTO/iTOsi,nr2 -to/^ -TOs^) . (15) 

As it is known, the separable force in Eq.® is expressed 
in the center of mass frame, the two-particle wave func- 
tion has to be transformed into the same frame. So we 
use the fact that the product of two oscillator functions 
in the coordinates of the two particles can be expanded 
in terms of products of oscillator functions in the relative 
and in the center of mass coordinates. This fact holds 
for the one-dimensional oscillators in z-direction 



as well as for the two-dimensional oscillators in r±- 
direction 



/ J / V ly pivipiLpriLp 

\plVlp ILplILp 



NpMp rip Trip 

Using the multinomial coefficients defined as 



Till 7712 ■■■ 'rui^ J mi!m2! . . . m^,! 



(16) 



with n = TOi + TO2 



TOjy for v — 2 (binomial coeffi- 



cients, which is expressed as and = 4. the one- 

\mi J 

and two-dimension Talmi-Moshinski transform brackets 
can be written as 1331: 
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M 



1 



(17) 



'"'^ V22^p+2"p+|Afpl+l™pl V (^p)!(^P + + |mp|)! 

E E E E(-)'^"^"* 



A^p A-/p rip mp 



Q,fl,S=0 T=0 g,r,s=0 t=0 



Np \ Mp \ Up 

Np~Q- R- S Q R S ) \ T J \np~q-r-s q r s 



r 



(18) 



Therefore the matrix element of the separable force in axially deformed oscillator basis can be written as: 

(12|y|l'2') = {7iz^nr^mi^,n^^nr^^mi^\V\nz^,nr^,mi^,,nz^,nr^,mi^,) = -G W^^^" ' Wy^!^" , (19) 



where 



with 



12 



(20) 



+00 



V"' - M 

"12 — -"^AT^ 



dx 



(21) 



^1? = aC7^:?7= 4 €M)e'^^-dx 



with the definition of 



%p = + Tlr^ + |m;j | - Np 



(22) 



and az — a/b^, ap — a/b±. Thus we find that the pair- 
ing matrix elements for the separable pairing interactions 
used in the RHB equation can be evaluated by a sum of 
separable terms in Eq. Using this expressions the 

pairing field A in Eq. ([5]) has the form 



ATp 



with 



Pn^N^ " I E ^12'^" '^12 



(23) 



(24) 



The results of the RHBZ model will depend on the choice 
of the effective RMF Lagrangian in the ph-channel, as 



well as on the treatment of pairin g co rrelations. In this 
work the effective interaction NL3 [33| is adopted for the 
RMF Lagrangian and in the pairing channel we use the 
separable form of the pairing force in Eq. ([8]) adjusted to 
the pairing part of the Gogny DIS force in Ref. [2l[. 



III. STUDY OF CONVERGENCE 

In the following investigations we solve the RHBZ 
equation (Eq. [3]) with this separable pairing force. As 
in Ref. [34] the Dirac spinors are expanded in an axially 
deformed oscillator basis with Np — 20 major oscilla- 
tor shells. As we see from the Eq. the separable 
pairing interaction is not fully separable in the axially 
deformed harmonic oscillator basis. We have a sum over 
the quantum number and Np characterizing the major 
shells of the deformed harmonic oscillator in the center 
of mass coordinate. For the self-consistent calculation 
of the RHBZ equation, the matrix elements of the V^^ 
and are calculated and stored in memory before the 
iteration. Therefore the time spend on the later calcu- 
lations of the pairing matrix elements is really negligible 
as compared to the total time, while it takes a large per- 
centage of computing time for the RHBZ model with the 
fully Gogny force in the pairing channel. Due to this big 
advantage much computer time can be saved in RHBZ 
calculations for axially deformed nuclei. 

For Np = 20 the maximum number of N^ and Np for 
the expansion of the pairing matrix elements in Eq. (|19p 
is N^ — 40 and iV^ — 20 respectively, which means a 
large number of 8000 separable terms in the pairing chan- 
nel. As we mentioned in last section, therefore a large 
memory space is required for this series. To reduce the 
storage we study the convergence of the expansion with 
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FIG. 1: (Color online) Neutron and Proton pairing energies 
of ^®'*Er obtained with different numbers of separable terms 
iV^^ and N° in Eq. JTIJ. 



the number of separable terms and and take the 
nucleus ^^''Er as example. For the self-consistent solu- 
tion of the RHBZ equations with the full expansion of 
the separable pairing matrix elements the deformation 
of the ground state of ^^"^Er is found to be prolate with 
P = 0.325. Varying the numbers and Np of the ex- 
pansion we plot in Fig. [T] the neutron and proton pairing 
energies of ^^"^Er with respect to the values of and 
Np. In the main panel of Figs, [ija) and [TJb), the con- 
tour lines correspond to lines of constant neutron and 
proton pairing energy. The value of these quantities in- 
creases with darkness. It is clearly seen that the neutron 
pairing energy converges for iV° = 14 and Np = 5. The 
top and right panels of Fig. [T{a) show how lihe neutron 
pairing energy converges as a function of N^ and Np for 
fixed values of Np = 5 and iV° = 14, respectively. The 
same is shown in Fig. [ijb) for the proton pairing energy. 
Due to the prolate shape of ^^"^Er the convergence in N^ 
is slower than that in Np. 

in Fig. IHa) the potential energy surface (PES) of the 
heavy nucleus ^"^^Pu is plotted as a function of the de- 
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FIG. 2; (Color online) The binding energy (a) and the pairing 
energy (b) as a function of the deformation for the nucleus 
^*°Pu obtained with different numbers of separable terms A^^ 
and 7V° in Eq. dHJ. 



formation parameter p. The full pairing matrix elements 
Eq. dUl) with iV° = 40 and 7V° = 20 are adopted in these 
calculations. The corresponding pairing energies of pro- 
tons and neutrons are given in Fig. [D^b) . It is found that 
the ground state of ^''"Pu has a normal prolate defor- 
mation (NpD) at /3npD — 0.278. Further two minima in 
the PES are observed with a normal oblate deformation 
(NoD) at /3noD = —0.277, and a super prolate deforma- 
tion (SD) at PsD = 0.865. As discussed above, a reduced 
number of separable terms with Np = 5 and N^ = 14 
in Eq. p9|) is large enough to obtain convergence for the 
pairing matrix elements at the normal prolate deforma- 
tion. We also investigate the convergence with the num- 
ber N^ and Np of separable terms for the other two min- 
ima in the PES. We find that for the super deformation 
a larger value of = 18 in the direction of the sym- 
metry axis is needed to obtain convergence for the global 
properties while in the perpendicular direction Np = 5 
is already large enough. On the other side, for the nor- 
mal oblate deformation, N^ = 10 is large enough, and we 
need at least Np = 7 to reach the full convergence. The 
details of these results are given in the Tab. [D 

In practical applications for normal deformed prolate 
nuclei it turns out that the expansion of the pairing ma- 
trix elements in Eq. (|19p can be restricted to finite values 
N;, < N^ = 14 and Np < N^ = 5 for obtain sufficient 
accuracy. For specific cases, where higher precision is 
required or for very large deformations, e.g. for superde- 
formed configurations larger values for N^ or Np are re- 
quired and convergence has to be checked. The matrix 



elements V-^" and V-^^' are calculated and stored before 
starting the iterations. As compared to the calculations 
with the full expansion of the separable force in the pair- 
ing channel this corresponds to a considerable reduction 
in the memory and computing time. 
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= -0.277 


E [MeV] 


En [MeV] 


Ep [MeV] 




full 




-1807.984 


-14.025 


-9.639 




10, Np = 


5 


-1807.123 


-8.870 


-8.718 




10, iVp = 


6 


-1807.819 


-13.022 


-9.609 




10, iVp = 


7 


-1807.972 


-13.975 


-9.606 


(3sD = 0.865 


E [MeV] 


En [MeV] 


Ep [MeV] 




full 




-1813.181 


-6.787 


-7.515 




14, iVp = 


5 


-1812.896 


-5.132 


-6.462 




16, iVp = 


5 


-1813.066 


-5.906 


-7.163 




18, iVp = 


5 


-1813.148 


-6.452 


-7.460 



TABLE I: The the total energy E and the pairing energies 
En and Ep for protons and neutrons for the nucleus ^^''Pu at 
normal and at super deformation for various values of A'^^ and 

n2. 
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the discrepancy of the binding energies obtained in these 
three models from the experimental data [s^. All the 
calculated binding energies are in good agreement with 
the experimental data. We observe that our results ob- 
tained from RHBZ-calculations with the new separable 
pairing interaction deviate from the experimental bind- 
ing energies by less than 0.2% and that they arc for many 
cases in slightly better agreement with experiment than 
the other models 

The deformations and shapes of nuclei play a crucial 
role in defining the properties such as nuclear sizes and 
isotopes shifts. They are strongly affected by the pairing 
correlations. In the Fig. |4] we show the quadrupole de- 
formation parameter /3 derived from Eq. ()13p for the Sm 
isotopes between A^=66 and A^=126. We find that the 
deformations obtained from RHBZ calculations with the 
new separable pairing force are very close to the results 
from non-relativistic HFB-calculations with the Gogny 
force DIS and in good agreement with the experimental 
data [lij. Although simple RMF + ECS theory provides 
a reasonable description for the binding energies, the de- 
formation parameters of neutron-rich nuclei calculated in 
this model deviate slightly from the experimental data in 
comparison with those obtained for RHBZ and HFB the- 
ory. 

Finally we show in Fig. [5] that there is excellent agree- 
ment between the calculations with the separable force 
presented in this paper and the original Gogny force in 
the pairing channel. The resulting deformation parame- 
ters shown in the upper panel are identical and for the 
pairing energies in the lower panel the differences are 
negligible. 
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FIG. 3: (Color online) Dependence of the binding energy (BE) 
and discrepancy between calculated binding energies and the 
available experimental data ABE on the number of neutrons 
for a chain of Sm(Z=62) isotopes. 



IV. SM-ISOTOPES 

We also perform RHBZ calculations with the separable 
pairing interaction for a chain of Sm (Z=62) isotopes in 
the rare-earth region. Several shape transitions are ex- 
pected along this isotope chain. The calculated ground 
state properties of Sm-isotopes, especially the deforma- 
tions are shown in Figs. [3] , [H and \E\ 

In the top panel of Fig. [3] we plot the total binding 
energy for Sm-isotopes as a function of neutron number 
for 66 < < 102. In comparison, we also show results 
obtained in two other theoretical models: (i) the RMF 
model with the parameter set NL3 [sst j in the ph-channel 
and the pairing correlations included by the Bardeen- 
Cooper-Schrieffer (BCS) formalism with constant pair- 
ing gaps obtained from the prescription of Ref. (36l | , and 
(ii) non-relativistic HFB calculations 37] with the Gogny 
force DIS In the lower panel of Fig. [3] we display 



V. SUMMARY 

We have presented first results obtained by axially 
symmetric relativistic Hartree-Bogoliubov calculations 
using the new separable pairing force introduced in 
Ref. This separable force is translational invari- 

ant and has finite range. It contains two parameters 
which are adjusted to reproduce the bell shape curve of 
the pairing gap at the Fermi surface obtained from the 
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FIG. 4: (Color online) Dependence of the deformation /3 on 
the number of neutrons for the chain of Sm isotopes between 
A=128 and A=188. 
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FIG. 5: (Color online) Comparison between RHB calculations 
based on NL3 with the original Gogny force DIS (stars) and 
with the new separable force adjusted to Gogny DIS (empty 
circles). Upper panel: Dependence of the deformation /3 for 
a chain of Sm isotopes between A=134 and A=154. Lower 
panel: the corresponding pairing energies 



of the pairing matrix elements converges relatively well 
and that an appropriate truncation provides an excellent 
approximation. This allows a considerable reduction of 
computing recourses such as time and memory in prac- 
tical applications of the RHBZ theory. In particular we 
study the ground state properties of a chain of well de- 
formed Sm-isotopes within this model. We find excellent 
agreement of our results with those obtained by using the 
non-relativistic HFB theory with the Gogny force DIS, 
the RMF -I- BCS theory based on the constant pairing 
gap approximation, and with the available experimental 
data. 

Results obtained by RHB-theory with the original 
Gogny force DIS and with the separable force derived 
from it are basically identical. Therefore we can con- 
clude that this simple pairing interaction can be applied 
in future applications of the RHBZ approach in nuclei far 
from stability instead of the complicated Gogny force. 
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